Inverse Energy Cascade in Forced 2D Quantum Turbulence 
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We demonstrate an inverse energy cascade in a minimal model of forced 2D quantum vortex turbulence. We 
simulate the Gross-Pitaevskii equation for a moving superfiuid subject to forcing by a stationary grid of obstacle 
potentials, and damping by a stationary thermal cloud. The forcing injects large amounts of vortex energy 
into the system at the scale of a few healing lengths. A regime of forcing and damping is identified where 
vortex energy is efficiently transported to large length scales via an inverse energy cascade associated with the 
growth of clusters of same-circulation vortices, a Kolmogorov scaling law in the kinetic energy spectrum over a 
substantial inertial range, and spectral condensation of kinetic energy at the scale of the system size. Our results 
provide clear evidence that the inverse energy cascade phenomenon, previously observed in a diverse range of 
classical systems, can also occur in quantum fluids. 



Turbulent fluid flows are characterized by the conserved 
transfer of kinetic energy across length scales (lj]. In clas- 
sical 3D turbulence, large-scale forcing leads to a cascade of 
energy and vorticity towards smaller scales 0]. In stark con- 
trast, small-scale forcing in 2D fluids may lead to an inverse 
energy cascade (IEC) Q3L |4£| , associated with energy flux to- 
wards larger scales, and the emergence of macroscopic ro- 
tating structures from the turbulent flow. The IEC has been 
widely studied in classical fluids J2, Si and observed over a 
large range of scales, from soap films S 01 to planetary at- 
mospheres fH). The IEC in classical fluids can be studied via 
a point- vortex model in which the growth of large-scale flows 
corresponds to the clustering of point-vortices with the same 
circulation JilQii. This model provides a link between classi- 
cal and quantum turbulence (QT), and suggests the intriguing 
possibility that the IEC may be observable in quantum fluids. 
The possibility to realize effectively 2D superfiuid flows has 
motivated the study of 2DQT in atomic Bose-Einstein con- 
densates (BECs) 1 11-22]]. and such flows have been experi- 
mentally demonstrated 123, 24]. Although intermittent few- 
vortex clustering has been observed in the breakdown of su- 
perfiuid flow around an obstacle JlH. 231. establishment of a 
steady, large-scale IEC in a BEC is hindered by the compress- 
ibility of the superfiuid, which allows vortex-antivortex anni- 
hilation lfl4lll5ll . and a direct cascade of acoustic energy via 
weak-wave turbulence |l2 ] . 

In this Letter we demonstrate an IEC of quantum vortices 
in numerical simulations of a forced, damped BEC. We show 
that the interplay between the effects of forcing and damp- 
ing leads to a regime where a clear IEC occurs. We analyze 
the dynamics of the 2DQT via the incompressible kinetic en- 
ergy (IKE) spectrum 115l 12511 associated with vortices, and the 
statistics of vortex clustering as identified by a cluster-finding 
algorithm. We observe three clear signatures of an IEC: (1) 
a sustained increase in the charge and spatial scale of vortex 
clusters; (2) a Kolmogorov kr 5 ^ IKE spectrum 12511 spanning 
an inertial range of up to 1.5 decades; and (3) spectral conden- 
sation lE^ - a macroscopic accumulation of IKE at the largest 
length scales of the system. 

Our system constitutes a BEC analog of forced 2D grid tur- 
bulence, a scenario that has been used to demonstrate the IEC 



in classical fluids |0]. We consider a superfiuid moving with 
initial speed v that is stirred by a stationary, regularly-spaced 
grid of four narrow Gaussian potentials. We model the dy- 
namics of the condensate, described by wavefunction ij/(r, t), 
using the damped Gross-Pitaevskii equation (dGPE) 12711 



#(r, t) 
dt 



(1) 



The operator _£ generates the zero-temperature GPE evolu- 
tion, and for a quasi-2D BEC subject to tight harmonic con- 
finement in the z-direction (with harmonic oscillator length L) 
can be written as 



-C = + V(x,y) + g 2 m,y, t)\\ 

2m 



(2) 



where V(x,y) is the grid potential, and g2 — 2 " [ /2jTn 2 a s /ml z , 
for atoms of mass m interacting with s-wave scattering length 
a s . Dissipative collisions between condensate atoms and a sta- 
tionary thermal cloud at chemical potential p. are described by 
the dimensionless damping rate y. Imposing periodic bound- 
ary conditions, our system provides a model of the obstacle- 
induced spin-down of apersistent current immersed in a sta- 
tionary thermal cloud 128 1. 

We work in energy, length and time units of chemical po- 
tential /i, healing length £ = H/mc, and g/c respectively, where 
c = ^fjTfm is the speed of sound in the homogeneous system. 
We numerically solve the dGPE in a periodic box of length 
L = 512£, for time t = 800^/c (at which time the turbulent 
wake reaches x = L/2). To eliminate the transient effects 
of abruptly raising the grid potential we begin from a ground 
state of the dGPE where the grid is co-moving with the flow, 
and stop the grid motion at t = 0. The grid potentials are 
then located at x — -L/2 + 8£ and equally spaced by L/4 
in the y direction. We add a small amount of initial noise to 
break the symmetry, and collect statistical information about 
the flow by computing multiple trajectories. Typical vortex 
configurations of the system at time t = 700^/c are shown in 
Fig-Q] To observe turbulent flow, we choose sufficiently high 
v « 0.822c such that the grid nucleates vortices in a chaotic 
manner fl^l . rather than periodically l29ll . To optimize vortex 
injection, we choose grid potential height Vo = lOOyU and 1 /e 2 
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FIG. 1. (Color online) Vortex configurations in a rightward-flowing superfluid stirred by a stationary grid of obstacle potentials, located at the 
left. Dots indicate the locations of positive (counter-clockwise circulation) and negative (clockwise circulation) vortices at time t = 700f/c. 
The vortices have been sorted into dipoles, free vortices, and clusters (see legend) using the cluster-finding algorithm (see text). Vortices in 
dipoles at the smallest length scales are not resolvable on the scale shown, and appear as isolated green dots. Streamlines (grey lines) give a 
visualization of the vortex-only velocity field by constructing the velocity field of an identical configuration of point vortices. 



radius wq = V8£- Each obstacle potential injects streams of 
opposite-signed vortices separated by distance ~ 2wq. As the 
linear momentum of a dipole is proportional to the vortex sep- 
aration, these streams carry negative momentum, arresting the 
superflow. 

The dimensionless damping rate y influences the vortex dy- 
namics, and can be controlled via the temperature. In decay- 
ing 2DQT, damping assists dipole recombination, and thus 
transfers IKE to small length scales via a direct energy cas- 
cade 111 511 . In our forced system damping is vital for generating 
an IEC. At low damping (Fig.[T| upper right panel), the vortex 
dynamics are dominated by mixing of the streams of opposite- 
sign vortices nucleated at each obstacle potential, leading to 
dipole annihilation. In contrast, high damping (Fig. [T] lower 
right panel) acts to reduce the linear momentum of the con- 
densate by rapidly polarizing (spatially separating) the vor- 
tex streams. This suppresses dipole annihilation until streams 
from adjacent obstacles merge. Intermediate damping (Fig.[TJ 
center panel) leads instead to the formation of small clusters 
of like-sign vortices. The subsequent merging and expansion 
of these clusters, ultimately leading to macroscopic dipoles 
approaching the scale of the system, is the spatial signature of 
an IEC in 2DQT JH. 

To quantify the vortex clustering seen in Fig. [TJ we have 
developed a cluster-finding algorithm that identifies the vor- 
tex clusters within a given configuration. This yields more 
detailed spatio-temporal statistical information regarding vor- 
tex clustering than scalar measures 12211 . The algorithm con- 
sists of two rules: (1) Opposite-sign vortices that are mutual 
nearest-neighbors constitute a dipole and are removed from 



the algorithm's consideration. (2) Same-sign vortices that are 
closer to each other than either is to an opposite-sign vortex 
are placed in the same cluster. Rule 1 (2) is applied recur- 
sively to all vortices still under consideration until no more 
vortices can be added to dipoles (clusters). Rule 1 is applied 
first, removing the short-range velocity fields of dipoles. The 
algorithm yields a decomposition of the configuration into 
clusters, dipoles, and free vortices that corresponds extremely 
well to the streamlines of the point-vortex velocity field, as 
seen in Fig. Q] 

In Fig.f2]we show the mean charge (number of vortices), k c , 
and radius (average distance of vortices from the cluster cen- 
ter), r c , of all clusters within a sample box of size (64^,5 12£). 
For each of the four values of y, we combine statistics from 
eight trajectories. In Fig. [2 a) we analyze the configurations at 
time t = 700£/c by centering the sampling box at (d - L/2,0), 
where d is the downstream distance from the grid. In Fig.|2jb) 
we analyze the configurations as a function of time, centering 
the sample box at (vf - 32^,0) [inset (vf - 96£,0)] such that 
it is co-moving with the flow. The results of these two analy- 
sis methods are in reasonable quantitative agreement, indicat- 
ing that the frozen-turbulence approximation (equivalence of 
d and vf) is satisfied 0]. The strong increases in k c and r c over 
large ranges of d and t in the case y = 0.009 are a clear sig- 
nature of an IEC. For lower values of the damping (y = 0.003 
and y = 0.0009) k c remains generally constant while r c in- 
creases slowly, suggesting the IEC is weak if present. For 
higher damping (y = 0.03) an initial rapid increase in k c , r c is 
followed by a similarly rapid decrease in both quantities; this 
is consistent with a strong IEC immediately behind the grid, 
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FIG. 2. (Color online) Spatio-temporal statistics of vortex clustering. 
For four values of y we show the mean charge, k c , and radius, r c , of 
clusters identified by the cluster-finding algorithm with cluster center 
located in a box of dimensions (64£,512<f), averaged over eight tra- 
jectories. Both quantities are shown as a function of: (a) downstream 
distance from the grid, d, using a box centered at (d - L/2,0) at time 
t = 700^/c; (b) time, t, using a box centered at (vf — 32£,0) and hence 
co-moving with the flow. Insets show: (a) distribution of positive 
(red plus) and negative (blue circle) clusters for the case y = 0.009; 
(b) k c and r c for a box centered at (vf - 96<f,0). By definition the 
smallest clusters consist of 2 vortices. 



which is destroyed by dipole recombination when the streams 
meet (Fig. [T}. 

The IKE spectrum of our system (Fig. [3} can be compared 
to the spectrum for an ideal model of 2DQT given in Ref . 121 ] . 
This model consists of quantum vortices with compressible 
cores; incompressible energy is efficiently transported from 
small-scale forcing towards large scales via vortex dynam- 
ics, but coupling of vortex dynamics to the sound field is ne- 
glected. In the ultraviolet (UV) region (k » the IKE 
spectrum is 



E\j{k) = CW)- 3 = A 2 M2£ 3 (^)-\ 



(3) 



where A is a dimensionless constant related to the structure 
of a compressible quantum vortex core, N is the total vortex 
number, Q = 2nH 2 no/m^ 2 is a constant with dimensions of 
enstrophy, and no = nlgi is the homogenous system density. 
The k~ 3 form of Eq. © is a universal feature of compressible 
2DQT arising from the finite size of a vortex core, and is not 
associated with a direct enstrophy cascade. At the wave num- 
ber k = there is a cross-over to point-vortex behavior in 
the infrared (IR) region (L _1 < k < k). Assuming (a) forcing 




FIG. 3. (Color online) Incompressible kinetic energy spectra at 
t = 550f/c for damping parameters (a) y = 0.0009 (b) y = 0.003 (c) 
y = 0.009 (circles), averaged over eight trajectories. Lines show the- 
oretical predictions for /T 3 [Eq. {3]l] and Kolmogorov k~ 5/3 [Eq. ©J 
spectra (see legend) in ideal 2DQT, with no fitted parameters. The 
Kolmogorov spectrum is shown using both the clustered vortex num- 
ber Nc (constant C*), and the total vortex number N (constant C) for 
comparison. 



wavenumberfc/- * k, and (b) all N vortices participate in clus- 
tering, the IR-region spectrum takes the Kolmogorov form 

E l c (k) = CW)- 5/3 (4) 

over an inertial range. The extent of the inertial range is linked 
to the scale range of vortex clusters, and it has been shown that 
appropriately distributed vortices in clusters of charge k c — 5 
are sufficient to generate a decade of inertial range. In our sim- 
ulations fey is determined by wq: treating the injected vortices 
as streams of dipoles of width // = 2 wo, forcing occurs near 
kf = 2n/lf, corresponding to the peak in the IKE spectrum 
of such a vortex dipole 12111 . Our choice wo = V&r thus cor- 
responds to kf = 1.1 Ik, satisfying assumption (a). Since the 
number of vortices participating in clustering, N c , is generally 
less than N, assumption (b) is not always satisfied. However, 
in this case the replacements Af — > N c and C — > C* = CN C /N 
in Eq. © give a better description of the inertial range, while 
retaining continuity with the UV region at k « k. 
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FIG. 4. (Color online) Fraction of the total kinetic energy, E tot , stored 
as incompressible kinetic energy at scale L (Eq, see text), as a func- 
tion of time. Inset: Vortex clustered fraction as a function of time. 



For all values of damping the UV spectrum agrees with the 
2DQT prediction [Eq. ©]. For y = 0.0009 [Fig. Ha)] the 
IR spectrum does not exhibit a clear k~ 5/3 scaling and lacks 
temporal stability. A buildup of energy around k is evident, 
suggesting a possible bottleneck ll30ll . associated with inef- 
ficient energy transport away from forcing wavenumber kj. 
For y = 0.003 [Fig. Ob)] the IR spectrum closely follows 
the modified 2DQT prediction [Eq. (0|i with N — > N c ] over 
an inertial range of around a decade of k, is steady for times 
300^/c < t < 600^/c, and has a weaker bottleneck around kf. 
This is consistent with increased vortex clustering (Fig. [TJ, 
and the onset of an IEC. For y = 0.009 [Fig. Etc)], the IR 
spectrum approximates the ideal 2DQT form [Eq. ©J over a 
wide inertial range (~ 1.5 decades) and is completely stable 
after t ss 400£/c. This is consistent with strong vortex clus- 
tering (Fig. [TJ, and a fully developed IEC. For strong damp- 
ing (y = 0.03) the spectrum remains qualitatively similar, al- 
though agreement with the ideal 2DQT IKE spectrum predic- 
tion is degraded by the abrupt surge in dipole annihilation oc- 
curring when vortex streams from adjacent obstacle potentials 
meet (Fig. U}. 

Flux of IKE to scales approaching the system size leads 
to an effect known as spectral condensation, associated with 
the emergence of macroscopic rotating structures |26|. Pre- 
viously, IKE fluxes in compressible, forced 2DQT have been 
calculated by neglecting the compressible part of the kinetic 
energy, or introducing assumptions regarding its coupling to 
the IKE spectrum 11141 1 1 311 - Here we provide a clear demon- 
stration of a negative IKE flux by considering the ratio of 
the IKE in mode k = 0, given by E = E\k = 0)Ak 
(where Ak = 2n/L), to the total kinetic energy, £ tot = 
(h 2 /2m) J d 2 r |V ± i/r(r)| 2 . We eliminate the contribution of the 
initial superflow to the IKE by computing the energies in the 
frame co-moving with the condensate at t = (Fig. |4j. For 



weak damping (y = 0.0009 and y = 0.003) the ratio Eo/E lot 
remains small. However, for higher damping (y = 0.009 and 
y = 0.03) the ratio E^/E tot initially grows rapidly, leading 
to significant spectral condensation. Merging of the vortex 
streams at later times (Fig.[TJ slows this growth (most notably 
for y = 0.03) and leads to a decline in the clustered fraction in 
all cases except for y = 0.009, where it remains consistently 
high (Fig. H] inset). Since the fractional occupation of the 
ground-state remains below 10~ 4 , this spectral condensation is 
not due to Bose-Einstein condensation into the ground-state, 
but rather emerges from macroscopic vortex clustering. 

We have demonstrated an inverse energy cascade in a quan- 
tum analog of classical 2D grid turbulence 0]. We solve the 
Gross-Pitaevskii equation describing a moving superfluid ar- 
rested by a grid of obstacle potentials, and damped by a sta- 
tionary thermal cloud. An inverse energy cascade is identi- 
fied via the growth and expansion of vortex clusters, a Kol- 
mogorov k~ 5 ^ incompressible kinetic energy spectrum, and 
a flux of incompressible kinetic energy to large scales. The 
demonstration of an inverse energy cascade of quantum vor- 
tices in a minimal model establishes a new link between clas- 
sical and quantum turbulence, and our methods of analysis 
form a foundation for future studies of quantum vortex tur- 
bulence in two dimensional superfluids. Our findings open 
new directions in the study of forced and decaying superfluid 
turbulence [16, 29| and provide a dynamical realization of the 
negative-temperature states for point vortices originally envis- 
aged by Onsager rf3~ill . 
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(UOO004). 
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